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Abstract 

The Fisher-Shannon information and a statistical measure of complexity are calcu- 
lated in the position and momentum spaces for the wave functions of the quantum 
isotropic harmonic oscillator. We show that these magnitudes are independent of 
the strength of the harmonic potential. Moreover, for each level of energy, it is found 
that these two indicators take their minimum values on the orbitals that correspond 
! to the classical (circular) orbits in the Bohr-like quantum image, just those with the 

• ' highest orbital angular momentum. 
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In recent years, the study of statistical magnitudes on quantum systems has 
taken an increasing interest [Tf2] . The tools that have been developed in in- 
formation and complexity theories, for instance, Fisher and Shannon informa- 
tions [3|^f5] . and statistical indicators of complexity [6|8] have been calculated 
for several systems under different approaches [9lll0|lll|ll2|ll3|ll4lll5j . The prob- 
ability densities characterizing the state of a quantum system are defined in 
the position and the momentum spaces [HUT]. From here, the calculation 
of all those statistical indicators can be performed with a low computational 
cost. 
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The dependence of these quantities on the quantum numbers of the system 
can reflect the hierarchical organization of that quantum system. Even for 
states with the same energy it is possible to have different values of these 
statistical magnitudes. Take, for instance, the i?-atom. It has been shown [TB] 
that for a given energy the minimum values of the Fisher- Shannon information 
and statistical complexity are reached for the highest allowed orbital angular 
momentum for that energy. This means that a variational process on these 
statistical measures can select just those orbitals that in the pre-quantum 
image are the Bohr-like orbits. 

Following this insight, it is our aim in the present work to analyze if the above 
described behavior of these statistical measures can be also found in the case 
of the isotropic quantum harmonic oscillator. 

Let us start by recalling the three-dimensional non-relativistic wave functions 
of this system when the potential energy is written as V(r) = A 2 r 2 /2, where 
A is a positive real constant expressing the potential strength. Atomic units 
are used through the text. The wave functions in position space (r = (r,Q), 
with r the radial distance and Q the solid angle) are: 

*nAm(r) = R n ,i(r) Y lim (p.) , (1) 

where R n> i(r) is the radial part and Y^ m (Vl) is the spherical harmonic of the 
quantum state determined by the quantum numbers (n, I, m). The radial part 
is expressed as [10J 



Rn,i(r) 



2,,!A ' +: "' J 1 ' 2 r l 4 +1/2 (Ar 2 ) , (2) 



T(n + I + 3/2) 



being L^(t) the associated Laguerre polynomials. The levels of energy are 
given by 

E n j = \(2n + I + 3/2) = \(e n j + 3/2), (3) 

where n = 0, 1, 2, • • • and / = 0, 1, 2, • • • . Let us observe that e n j = 2n + I. 
Thus, different pairs of (n, I) can give the same e n j, and then the same energy 
E n ,i- 

The wave functions in momentum space (p = (p, U), with p the momentum 
modulus and Cl the solid angle) are: 

*n,I,m(p) = RnM Y lA^) , (4) 

where the radial part R n j(p) is now given by the expression [10J 



Rn,l(p) 



T(n + I + 3/2) 



p l e-h L l n + ^(p*/\) . (5) 
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Taking the former expressions, the probability density in position and mo- 
mentum spaces, 

Pn,l,m;x(r) = | ^ n ,l,m(r) | 2 , 7n,l,m; \(p) = \ ^ n ,l,m{p) ? , (6) 

can be explicitly calculated. From these densities, the statistical complexity 
and the Fisher- Shannon information are computed. We find that these quanti- 
ties are independent of A, the potential strength. This is a non trivial property 
that is shown in the Appendix A. For this reason, we drop the A subindex in 
the densities from now and on. Also, for the sake of simplicity, the quantum 
numbers (n, I, m) are omitted in the notation. 

First, the measure of complexity C recently introduced by Lopez- Ruiz, Mancini 
and Calbet [6]T7t8] . the so-called LMC complexity, is defined as 

C = H-D , (7) 

where H represents the information content of the system and D is the distance 
from the actual state of the system to some prestablished reference state. 

For our purpose, we take a version used in Ref. [S] as quantifier of H. This is 
the simple exponential Shannon entropy, that in the position and momentum 
spaces takes the form, respectively, 

H r = e Sr , H p = e Sp . (8) 

S r and S p are the Shannon information entropies [I], 

S r = - J p(f) logp(f) dr, S p = - J 7(p) log7(j9) dp. (9) 

We keep for the disequilibrium the form originally introduced in Refs. [6f8] . 
that is, 

D r = J P 2 (r) df, D p = J 7 2 ($ dp, (10) 

In this manner, the final expressions for C in position and momentum spaces 
are: 

C T = H r ' D r , Cp = Hp ' Dp . ( 1 1 ) 

The form of the wave functions, due to the harmonic interaction, allows us to 
show in the Appendix A that these magnitudes, C T and C p , are the same. 

In Fig. 1, C r (or C p ) is plotted as function of the modulus of the third com- 
ponent m, —I < m < I, of the orbital angular momentum / for different I 
values with a fixed energy. That is, according to expression ([3]), the quantity 
e n j = 2n + I is constant in each figure. Fig. 1(a) shows C r for e n j = 15 and 
Fig. 1(b) shows C r for e n ,i = 30. In both figures, it can be observed that C r 
splits in different sets of discrete points. Each one of these sets is associated 
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to a different I value. It is worth to note that the set with the minimum values 
of C r corresponds just to the highest I, that is, I = 15 in Fig. 1(a) and I = 30 
in Fig. 1(b). 

Second, other types of statistical measures that maintain the product form of 
C can be defined. Let us take, for instance, the Fisher- Shannon information, 
P, that has been also applied in Refs. [TT1I15II18] for quantum systems. This 
quantity, in the position and momentum spaces, is given respectively by 

P r J r • I r , Pp Jp • Ip , (12) 

where the first factor 

J r = — e 2S ^ 3 , J p = — e 2S ^ 3 , (13) 
is a version of the exponential Shannon entropy [5J, and the second factor 



W**-, , S= /M dj r, (14) 



p(f) ' " J 7(p) 

is the so-called Fisher information measure [3] , that quantifies the narrowness 
of the probability density. Similarly to the behavior of C r and C p , we also 
show in the Appendix A that P T = P p . 

I r can be analytically obtained in both spaces (position and momentum). The 
results are |12|: 



J r = 4 (2n + / + 3/2 - \m\) X , (15) 
Jp = 4 (2n + I + 3/2 - \m\) A -1 . (16) 

Let us note that I r and I p depend on A, although the final result for P r and 
Pp are non A-dependent (see Appendix A). 

Fig. 2 shows the calculation of P as function of the modulus of the third 
component m for different pairs of (e n> j = 2n + l,l) values. In Fig. 2(a), P r (or 
P p ) is plotted for e n j = 15, and P r is plotted for e n ^ = 30 in Fig. 2(b). Here, 
P r also splits in different sets of discrete points, showing a behavior parallel 
to the above signaled for C (Fig. 1). Each one of these sets is related with a 
different I value, and the set with the minimum values of P r also corresponds 
just to the highest I, that is, I = 15 and I = 30, respectively. 

Let us finish this note with the conclusions. The statistical complexity and 
the Fisher-Shannon information have been shown to be independent of the 
potential strength, A. It is the specific forms, (IB]) . (flOl) and (fT3j) . (fT4|) . of the 
definitions of these two indicators that yield this striking property. This fact 
could be an indirect argument to justify the choice of these expressions. Then, 
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these magnitudes have been calculated. We have taken advantage of the exact 
knowledge of the wave functions. Concretely, we put in evidence that, for a 
fixed level of energy, let us say e Ui i — 2n + l, these statistical magnitudes take 
their minimum values for the highest allowed orbital angular momentum, I = 
e n j. It is worth to remember at this point that the radial part of this particular 
wave function, that describes the quantum system in the (n — 0, 1 — e n> /) 
orbital, has no nodes. This means that the spatial configuration of this state 
is, in some way, a spherical-like shell. In the Appendix B, the mean radius of 
this shell, < r > ni ; i?n , is found for the case (n — 0, 1 — e n j, m). This is: 

< r > n=0 ,i= e „ , Jtm =< r > n=0) i =e „ ii ~ y/ X- 1 (e n> i + 1) (l + e(e~j)) , (17) 

that tends, when e n j ^> 1, to the radius of the Nth energy level, = 
JX" 1 (N + 1), taking N = e n j in the Bohr-like picture of the harmonic oscil- 
lator (see Appendix B). 

As it was remarked in Ref . [18] , here we also obtain that the minimum values 
of the statistical measures calculated from the wave functions of the quantum 
isotropic harmonic oscillator select just those orbitals that in the pre-quantum 
image are the Bohr-like orbits. Therefore, we conclude that our intuition is 
enhanced when using these magnitudes to discern complexity at a quantum 
level. 



APPENDIX A: Invariance of C and P under rescaling transforma- 
tions 



Here, we show that the statistical complexities C r and C p are equal and in- 
dependent of the strength potential, A, for the case of the quantum isotropic 
harmonic oscillator. Also, the same behavior is displayed by P r and P p . 

For a fixed set of quantum numbers, (n,l,m), let us define the normalized 
probability density pit): 

M= nn™ 3 /2) *°-*K"W'^™ 2 - (18) 
From expressions (jTJ), (j2J) and (jSJ), it can be obtained that 

p A (f) = A 3/2 p^r), (19) 

where p\ is the normalized probability density of expression (Jf)]). Now, it is 
straightforward to find that 

H r {p x ) = A- 3 / 2 H(p) , (20) 
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and that 

D r (p x ) = A 3 / 2 £>(/>)• (21) 

Then, 

C r (p x )=C(p), (22) 
and the non A-dependence of C r is shown. 

To show that C r and C p are equal, let us note that, from expressions (TjJ, ([5]) 
and ([6]), the normalized probability density 7a (p) for the same set of quantum 
numbers (n, I, m) can be written as 

7A(p) = A- 3 / 2 p(A- 1 / 2 p). (23) 



Now, it is found that 

and that 

Then, 



H p { lx ) = A 3/2 H{p) , (24) 
D p { lx ) = A- 3 / 2 D{p). (25) 



C p ( 7 a) = C(p), (26) 
and the equality of C r and C p , and its non A-dependence are shown. 

Equivalently, from expressions ffl~2l . ([TBI . ffT5l) and ([16]) . it can be found that 
P r = P p , and that these magnitudes are also non A-dependent. 



APPENDIX B: Bohr- like orbits in the quantum isotropic harmonic 
oscillator 



Here, the mean radius of the orbital with the lowest complexity is calculated 
as function of the energy. Also, the radii of the orbits in the Bohr picture are 
obtained. 

The general expression of the mean radius of a state represented by the wave 
function ^> n) i )m is given by 



< r > n r m =< r > n r- 



711 



T(n + / + 3/2) A 1 / 2 Jo 



t l+1 e-* 



L l + 1/2 (t) 



dt. 



(27) 

For the case of minimum complexity (see Fig. 1 or 2), the state has the quan- 
tum numbers (n = 0, / = e n j). The last expression (1271) becomes: 



(eni + 1)! 



: e n ,l" 



r(e n , + 3/2)AV2 



(28) 
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that, in the limit e n> ; 3> 1, simplifies to expression (ITT)) : 

< r >ri= o iJ=eB>I>1 ~ ^A-i (e n>l + 1) (l + 9(e-J)) . (29) 

We proceed now to obtain the radius of an orbit in the Bohr-like image of 
the isotropic harmonic oscillator. Let us recall that this image establishes the 
quantization of the energy through the quantization of the classical orbital 
angular momentum. So, the energy E of a particle of mass m moving with 
velocity v on a circular orbit of radius r under the harmonic potential V(r) = 
m\ 2 r 2 /2 is: 

E = Kn\ 2 r 2 + ^mv 2 . (30) 
The circular orbit is maintained by the central force through the equation: 



mv 2 

= mX r. (31) 

r 

The angular momentum takes discrete values according to the condition: 

mvr = (N + l)h (N = 0,1,2,...). (32) 



Combining the last three equations (1301) -(1321). and taking atomic units, m 
H — 1, the radius of a Bohr-like orbit for this system is obtained 



r N 



^\-\N+l) (N = 0,1,2,...). (33) 



Let us observe that this expression coincides with the quantum mechanical 
radius given by expression (129]) when e Ui i = N for N ^> 1. 
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(a) (b) 

Fig. 1. Statistical complexity in position space, C r , and momentum space, C p , vs. 
I ml for different energy e ni /-values in the quantum isotropic harmonic oscillator for 
(a) e n ,l = 15 and (b) e n j = 30. Recall that C r = C p . All values are in atomic units. 
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Fig. 2. Fisher-Shannon information in position space, P r , and momentum space, P p , 
vs. \m\ for different energy 6 n? /-values in the quantum isotropic harmonic oscillator, 
for (a) e n< i = 15 and (b) e n j = 30. Recall that P r = P p . All values are in atomic 
units. 
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